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This paper is concerned with string cosmology and the dynamics of multiple scalar fields in 
potentials that can become negative, and their features as (Early) Dark Energy models. Our point 
of departure is the "String Axiverse" , a scenario that motivates the existence of cosmologically light 
axion fields as a generic consequence of string theory. We couple such an axion to its corresponding 
modulus. We give a detailed presentation of the rich cosmology of such a model, ranging from the 
setting of initial conditions on the fields during inflation, to the asymptotic future. We present some 
simplifying assumptions based on the flxing of the axion decay constant fa, and on the effective 
field theory when the modulus trajectory is adiabatic, and find the conditions under which these 
assumptions break down. As a by-product of our analysis, we find that relaxing the assumption of 
fixed fa leads to the appearance of a new meta-stable de-Sitter region for the modulus without the 
need for uplifting by an additional constant. A dynamical systems analysis reveals the existence 
of many fixed point attractors, repellers and saddle points, which we analyse in detail. We also 
provide geometric interpretations of the phase space. The fixed points can be used to bound the 
couplings in the model. A systematic scan of certain regions of parameter space reveals that the 
future evolution of the universe in this model can be rich, containing multiple epochs of accelerated 
expansion. 



I. INTRODUCTION 

String cosmology [iHl] has been undergoing a renais- 
sance over the last decade, due in part to our increased 
understanding of the landscape [S', , advances in mod- 
uli stabilisation and string model building ^-Q^ and the 
continued successes of precision cosmology [lOj. Theory 
and experiment have found fruitful harmony in two main 
areas: the early time accelerated expansion of the uni- 
verse and generation of primordial density perturbations 
during inflation fTTHTT] and in the study of the current 
epoch of cosmological acceleration [THl US] through the- 
ories of dark energy and quintessence [30]. Recently, the 
"String Axiverse" [21] has motivated the study of ultra- 
light axion fields that evolve on cosmological time scales 
P^l - BS] : one should ask whether the axion dynamics in 
this scenario can have other interesting or adverse effects 
on string cosmology and model building. 

Axions have been ubiquitous in theoretical physics 
since they were first proposed to solve the strong CP 
problem [2S]. Ever since this time, they have presented 
many problems and possibilities to cosmologists [SDHIS] 
(for reviews of axion cosmology, see [46l |47j). Indeed, 
axions are the major contender with Weakly Interacting 
Massive Particles (WIMPs) as a main constituent of the 
Dark Matter (DM). If both axions and WIMPs exist, it 
is natural that they should share the DM burden J48; . 

In this paper we systematically explore the parameter 
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space and dynamics of the coupled axion-modulus system 
originally proposed in |49| . which is a phenomenologi- 
cal extension of the axiverse to include effects on mod- 
uli. The energy scales of this model cannot be argued so 
model independently and elegantly to give rise to inter- 
esting cosmology as in the axiverse alone. However we 
take a phenomenological approach and use it to ask the 
questions: can cosmological axion dynamics lead to cos- 
mological modulus dynamics? What are these dynamics? 
For what parameter values does interesting phenomenol- 
ogy occur? It is the purpose of this paper to answer these 
questions. 



When dimensionally reduced and considered at low en- 
ergies, string theory furnishes us with extra degrees of 
freedom, in addition to the ACDM concordance model in- 
gredients of General Relativity, cold dark matter (CDM), 
and the standard model of particle physics. In fact, cos- 
mologically relevant axions and moduli are the generic 
prediction of string/M-theory 50 . These extra degrees 
of freedom can be viewed as a blessing or a curse. Ex- 
tra scalar fields with appropriately fine tuned potentials 
in the early universe are useful for inflationary model 
building, but they also lead to the cosmological moduli 
problem and must be properly stabilised. Scalar fields 
can also serve as dark matter, or dark energy. The sim- 
plest 6-parameter version of concordance ACDM cosmol- 
ogy may soon be observationally extended with the de- 
tection of neutrino mass and mass splittings [51] ■ There 
are also observational hints from small-scale CMB exper- 
iments |52j 53 that other relativistic degrees of freedom 
can already be seen in cosmology. The possibility that 
this may be string theory related, for example through 
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Symbol 


Meaning 


Reference Eq. 


fa 


axion decay constant 




1 




Aa 


axion potential energy scale 




1 




M 


energy scale of non-perturbative physics for axion potential 




2 




B,D 


modulus potential parameters 




4 




PA 


cosmological constant energy density 




4 




A 


value of the true vacuum energy 


( 


24 


1 




axion field 




3 




X 


modulus field 




4 




c 


axion-modulus coupling constant 




4 




7b 


baryotropic fluid equation of state 




6 






autonomous system variables 


( 


19 


1 


iVcfd 


number of e-folds from beginning of model evolution until end of fluid domination 




A/'ae 


number of periods of accelerated expansion 




fie 


early dark energy (EDE) density 


( 


42 


1 


M 


— n'^/fa (axion mass scale) 




4 




/? 






12 






— VG/P (critical value of (j> above which modulus is destabilised) 




13 






= VB/D 




11 




c 


= Bpa/D^ = t^z^/r'^ (vacuum constraint). C > (<)l/4 gives dS (AdS) 




23 





TABLE I: Symbols used in this paper. 



models of Early Dark Energy [MlES], is an exciting one. 
Independently of string theory, a multi-component dark 
sector is natural in, and may be strong evidence for, the 
top-down or anthropic view of cosmology [48] . 

The landscape of string theory vacua can also be seen 
as a blessing or a curse: it hampers the exact predictiv- 
ity of string theory, but it may yet explain the smallness 
of the cosmological constant by providing a high enough 
density of possible vacua "near" to a phenomenologically 
acceptable one [B], and through eternal inflation provide 
a mechanism for the scanning. In the model we will 
study, the overall value of the vacuum energy is a free 
parameter; it was shown in [49] that one can have this 
negative, so that the universe today is rolling towards 
collapse. Harlow et al [56, 57j have argued that the exis- 
tence of terminal vacua is necessary for the existence of 
a global arrow of time in eternal inflation, and studying 
their phenomenology therefore seems pertinent. In all 
cases what we gain from string cosmology is the poten- 
tial to ask deeper questions when models are embedded 
in a UV-complete theory. 

In our view, another bright side of the string land- 
scape comes from looking for the "why not?" features. 
The axiverse is one such generic feature. Along with the 
moduli of the landscape, we also get axions. We will 
save the details of this scenario for later, for now all we 
need say is that many of these axions should remain light, 
and thus cosmologically active. Through their effects on 
structure formation |23) it will be possible with next- 
generation cosmological observations to constrain the ex- 
istence of these fields as a component of the dark sector 



energy density at percent-level accuracy [26]. The axi- 
verse raises the possibility that string cosmology may be 
active at late times and that we may be able to observe 
it. This has been called the "Low- Energy Frontier of 
Particle Physics" [58] . If we can detect axions with high 
precision using cosmology, might we also detect changes 
in their evolution caused by the moduli? 

As the simplest models start to be constrained, we can 
begin to explore them more deeply. Naturally, one goes 
from assuming that just one scalar field is active at late 
times, to assuming that many are f59H61| . Already in 
inflation, thought of as embedded in string theory, the 
dominant paradigm is of an inflationary direction in a 
multi-dimensional field space (e.g. [B^HSl])- It should 
also be true that the light axions of the axiverse exist 
as a flat direction in some much larger field space of 
their brother axions and sister moduli. To look simply 
at many axion fields, and ignore their partner moduli is 
the most conservative option: the axiverse should be gen- 
eral enough to exist independently of the mechanism for 
moduli stabilisation. But we will show here that there 
are cases where one cannot ignore the effects of the mod- 
uli, and so we explore what these effects might be within 
a certain parameterisation. 

This paper is organised as follows: In Section ([ll| we 
recall the coupled axion-modulus system introduced in 
[32] and comment on how the initial conditions appropri- 
ate for acceptable late time phenomenology might arise. 
In Section ( III ) we analyse the basic features of cosmol- 



ogy in this model by looking at the scalar potential, be- 
fore presenting a systematic analysis of the dynamics us- 
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ing a dynamical systems approach in Section (IV I. We 
discuss our findings in Section (Tv]), and conclude in Sec- 



tion (VII 



Our detailed presentation of the rich cosmology of our 
model requires the use of many equations and symbols. 
As a guide to the reader, in Table ^ we summarize the 
key symbols and the equations where they are defined or 
first used. 



II. THE MODEL 



to some ultra-light, stable axions for the cosmologist to 
play with. 

We canonically normalise (j) = faO, and expand the 
potential around the minimum: 



C 



2 K 
J a 



(3) 



The Axiverse 

All types of string theory and M-theory contain mul- 
tiple axion fields [551 ES] ■ These axions arise when anti- 
symmetric tensor fields are compactified on closed cycles; 
the axion is the Kaluza-Klein zero mode, and appears in 
the gauge kinetic function. Axions then acquire a po- 
tential when non-perturbative physics is turned on on 
the cycle, for example from wrapped D-branes, or from 
world sheet or gauge theory instantons. The axion is 
a pseudo-Nambu-Goldstone boson (PNGB) of a spon- 
taneously broken global symmetry. Many extensions of 
the standard model of particle physics also contain more 
generic PNGBs [67|. 

The low energy four dimensional Lagrangian for an 
axion, 9, with periodic potential U {6) is: 

c^^idof-Atuie). (1) 

The two scales in this Largrangian, the decay constant 
fa, and potential energy scale A^, both depend upon the 
action, 5*, of the non-perturbative physics on the corre- 
sponding cycle in the following way: 

S' ' 

(2) 

Here Mpi is the (reduced) Planck mass and fj, sets the 
scale of non-perturbative physics, for example the QCD 
scale or, in string theory, the geometric mean of the 
supersymmetry (SUSY) breaking scale, Mgugy and the 
Planck scale, ^^/Mp^Mrttsv-^PT]. The action in turn de- 
pends on the size of the cycle, and herein we find the 
axiverse mechanism for light axions. Although fa should 
be fixed at some high scale fa ~ 10^^ GeV |66j, small 
variations in the sizes of cycles and the exponential sen- 
sitivity of Aa on S means that the axion mass should 
distribute roughly evenly on a logarithmic scale, leading 



^ In |49| this mean was misquoted. A forthcoming erratum will 
correct this. The argument and scales in that work should follow 
from a higher scale of SUSY breaking 



The Axiverse Is More Than Just Axions 

As with everything in string theory, our low energy 
"constants", such as nia, are not really constant at all, 
but depend upon moduli. In this case the modulus of 
interest is that controlling the area of the cycle giving 
us the axion. S depends on this area, and so we can 
choose to identify S = Cx for some modulus field x, and 
coupling C. Eq. ^ then implies that axions and moduli 
are coupled. 

This modulus must be stabilised non-perturbatively 
|68) . Typical potentials are sums of exponentials (see, 
for example, I69j ). These considerations led, in 

[H] , to the following potential being studied for a coupled 
axion- modulus system^: 

F(0, x) = Be-^^^ ~ De-^^ + ^^M^e'^^cf" + pa , (4) 

where pA is the cosmological constant, added arbitrarily 
in this model so that whatever value the potential takes 
today can be made consistent with observations. We also 
have — n'^/fa- The total Lagrangian is then of the 
form: 

C=lid^r + lidxf-V{cl^,x). (5) 

There is one important caveat to this picture: by tak- 
ing fa fixed we are implicitly assuming small modulus 
variations. This greatly simplifies our system, since if 
we allowed for the variation of faix) this would change 
the canonical normalisation of the axion kinetic terms, 
and introduce kinetic mixing between the (f) and x fields. 
This effect could introduce new phenomenology in ex- 
treme trajectories with large Ax/Xj for example the pos- 
sibility of chaotic behaviour, but we defer study of this to 



^ (i) The general form of the potential we study fits into the class 
of models of "Generalized Assisted Inflation" I70H72| . though 
we mainly emphasise its use for quintessence, rather than in- 
flationary, purposes, (ii) At large values of Xi loop effects will 
eventually cause the potential to rise again, as in [9l l73| . We do 
not consider such contributions. Our conclusions only depend on 
having a sufficiently long, flat region of the potential before these 
effects kick in, but we do not compute the scales of parameters 
necessary for this. 
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a future work^. We make some comments on the effect 
on the potential in Section ( [v[ ). 

We assume that the axion and modulus fields evolve 
in a spatially flat Friedmann-Lemaitre- Robertson- Walker 
(FLRW) background containing the bare cosmological 
constant, pA, and a single fluid with baryotropic equa- 
tion of state Pb = (7b — l)Pbj where 7b is a constant, 
< 7b < 2. For radiation 7b = 4/3 or for dust CDM 
7b = 1. In standard cosmic time the evolution equations 
are: 

• dV 

dV 

Pb = -3i?pb7b , 

H = -^[Pb7b + </>'+x']- (6) 
These are subject to the Friedmann constraint 

= Pb + pA + + Px^ (7) 

where we have assumed that we can define the following 
distinct densities and pressures: 

P4> - + u{4>,x). 

Px = Ix' + Vb{x)-Vd{x), 
Px = \x^-Vb{x) + Vd{x)- (8) 
by splitting the potential as: 

V{c^,x)^Vb{x)-Vd (x) + C/(0, x) + pa , (9) 

where 

Vb{x) = Se-2f^^ 
Vd{x) = De-^^, 

U{C^,X) = \e-''^M^4,\ (10) 

We note that the split in these cases between dark mat- 
ter and dark energy is somewhat arbitrary |74| . When 
scalar fields begin oscillating, they redshift and cluster 
as dark matter, with an individual equation of state 
W(jix ~^ Oj but the modulus term also has a vev, behaving 
as a cosmological constant. Slow roll, and the inclusion 
of Pa muddy the waters further. Also, it is completely 
arbitrary in the presence of coupling to place U{4>,x) i^i 
P4>- 



^ We thank John March-RusscU for pointing this fact out to us. 



In the axiverse as presented in all moduli were as- 
sumed to be absolutely stabilised. The different sizes at 
which they were stabilised led to the different masses for 
the axions. Some were stabilised at larger values than 
others in order to make some axions light, but the differ- 
ences are not hierarchical. Assuming absolute stability 
implied that the moduli were heavy, and lived in their 
global minimum. The axiverse has been concretely re- 
alised in the moduli stabilisation scheme of [32]. In the 
scheme of 0, moduli were stabilised at hierarchically dif- 
ferent values, which allowed some moduli to remain very 
light. In this picture the axion phenomenology is not 
considered, and they are set to their vacuum values at 
zero. This is perfectly well justified even for light axions 
if all we are concerned with is the existence of a stable 
minimum for the moduli and in calculating their masses 
at this minimum, but if the axion evolution is our focus, 
then their possible effects on the moduli cannot be ig- 
nored. In both [22j and [9j only a handful of fields were 
considered, not the hundreds motivated in the axiverse. 
Here we take inspiration from the success of these models, 
and apply to it the spirit of optimism of the axiverse to 
look for phenomenology in a larger arena of possibilities. 



Cosmic Overview 

Here we give an overview of a scenario that may lead to 
the realisation of the initial conditions appropriate to our 
model, and the picture of cosmic history that emerges. 

We will assume that the universe begins in an eter- 
nally inflating de-Sitter (dS) false vacuum. This vacuum 
decays via tunnelling and bubble nucleation |75| into the 
standard phase of slow-roll inflation required to generate 
the primordial power spectrum. It must also be assumed 
that initial conditions on the axion and modulus fields 
are laid down prior to inflation. After slow-roll inflation 
ends, the inflaton decays and reheats the universe. We 
will consider the toy model of a post-inflation universe 
consisting only of matter, radiation, and the axion and 
modulus field condensates contributing a dark sector en- 
ergy density as described in [49] . 

In [IH], initial conditions were such that the modulus 
began at a large value, Xi and the axion mass rria = 
M'^e~^^ was cosmologically hght, = 0(1 - IQ^°)Hq. 
This modulus initial condition was not at the local min- 
imum of the potential, (x)(0). The axion initial condi- 
tion, (j)i is set at the Peccei-Quinn (PQ) phase transition 
by spontaneous symmetry breaking [291 176] . There are 
two logical possibilities for modulus evolution: there is a 
local modulus minimum at 0^, or there is not. If there 
is a minimum, and the modulus is sufficiently heavy to 
overcome Hubble friction, it will, like the inflaton before 
it, roll to the local minimum, while the light axion frozen 
at prevents it from reaching the global minimum. If 
there is no local minimum, then the modulus will roll to 
yet larger values until stopped by Hubble friction, only 
decaying to a local minimum once one exists (see e.g. 
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[77])^. Under these conditions, although the modulus 
and axion masses at the global minimum could be large 
(e.g. string/Planck scale), interactions instead freeze the 
modulus either in its local minimum, or at x > Xi- 
both cases, the modulus must evolve with the axion. 

For these initial conditions to be possible the only re- 
quirement is that the PQ phase transition and the switch- 
ing on of the appropriate instanton effects, which create 
the axion condensate and form the coupled potential, 
happen before the modulus finds its minimum. Certainly, 
during slow-roll inflation there are scalar fields yet to find 
their minima: the inflaton itself is one such field, ft is 
not unreasonable to assume that there are other moduli 
present that also exist away from their minima. Indeed, 
this is the case in any model of multi-field inflation and is 
the string interpretation of any quintessence model. This 
is also the expectation for the post-inflation, pre-hot-big- 
bang phase in string cosmology, where the post-inflation 
universe is dominated by the yet-to-decay, matter-like 
moduh gum HO]. 

The requirement that inflation occurs after the PQ 
phase transition, and that the reheat temperature does 
not restore the PQ symmetry, is generic to almost all 
models with axions as it is required to avoid a cosmolog- 
ical abundance of disastrous domain walls and the like 
[76j . In addition, string axion models require a low en- 
ergy scale of inflation [21 , in part to avoid overproduc- 
tion of axion isocurvature perturbations. 

After the radiation and matter dominated phases end, 
the next stage in the evolution of the universe again has 
a number of possibilities depending on the axion and 
modulus flelds. The additional cosmological constant in 
the potential, p\, can be regarded as the usual left over 
contribution to the vacuum energy. It has contributions 
which reduce it over time after spontaneous symmetry 
breaking (e.g. at the electroweak transition), positive 
contributions from dS stabilised moduli and vacuum fluc- 
tuations of quantum flelds, and negative contributions 
from AdS stabilised moduli (see e.g. [6]). The value of 
this constant relative to the potential minimum in the 
axion and modulus flelds determines the fate of the uni- 
verse. 

If the axion and modulus flelds are heavy and their held 
values are small, such that they are oscillating about, 
or slowly rolling into, their global minimum during the 
present epoch, then the scenario will be much like any 
other quintessence or axion dark matter scenario. The 
value of the total cosmological constant in the bottom 
the potential must be small, and of the correct magni- 
tude to account for the observed accelerated expansion 
of the universe. However, if the flelds are light enough 
and their initial values large enough that they are on the 



plateau of the potential, then the phenomenology can be 
quite different. Here, the potential energy of the axion 
and modulus are a small contribution, and the current 
accelerated expansion will be driven almost entirely by 
PA, as was the case in [15] . 

However, as also described in [l^, the axion-modulus 
system is only quasi-stable: eventually axion oscillations 
will decay and the modulus will flnd the global minimum. 
Depending on the initial conditions and the parameters 
in the potential it is possible to arrange for an acceptable 
cosmology where the vacuum energy at the global min- 
imum is either positive, negative, or indeed zero. If the 
vacuum energy is negative then the decay of the modulus 
will trigger rapid cosmological collapse, rather than life 
in a stable AdS state [75J [511 [SI]- In the case where it 
is positive, then a scenario such as explored in [53] will 
ensue. 



III. AXION-MODULUS DYNAMICS AND THE 
COUPLED POTENTIAL 

Before we begin our detailed dynamical systems anal- 
ysis in Section ( IV ) , we aim to give here some basic intu- 



ition about the types of phenomena possible in a cosmol- 
ogy with coupled scalar fields and an arbitrary vacuum 
energy. In particular, we give examples of phenomena 
not explored previously in [49 . The examples use ar- 
bitrary values of the parameters and are for illustration 
only. 



Local Minima, the Adiabatic Trajectory, and an 
Effective Potential 



The form of the potential is such that there is just 
one minimum, when the axion is at zero. However, for 
light axions, most of cosmic history ^ is spent away from 
this global minimum. The form of the coupling between 
axion and modulus then means that the moduli, too, will 
live away from their global minimum and consequently 
moduli must evolve during the course of cosmic history. 

We find that the local minimum in the modulus direc- 
tion as a function of the axion background is: 



1 

"c 



In 



1 

2^ 



1 



where 



3M2 



(11) 



(12) 



* These considerations are basically a statement of the Brustein- 
Steinhardt problem |78) for this model. Related issues are dis- 
cussed in |79| . 



When viewed in redshift space, where we have cosmic coincidence 
for things occurring within 2 = 1 of us, which is of order billions 
of years. The coincidence problem is in fact much more of a 
problem forwards in time: why are we not Boltzmann brains in 
thermal de-Sitter space? 
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FIG. 1: The potential of Eq. Q near the global minimum, 
for arbitrary parameters. In red, the adiabatic trajectory: 
the modulus minimum as a function of (Eq. (111). When 
(j) takes large values (defined by Eq. (13l), the minimum at 
finite X is destroyed. In purple, the naive trajectory: the 
axion potential at fixed x ~ {x)(0). For a heavy modulus the 
adiabatic trajectory will be followed, which is shallower near 
the minimum than the naive trajectory. 



and oj^ = B/D. We plot this trajectory on the potential 
surface in Fig. ([T]). 

The existence of the global minimum at positive mod- 
ulus translates into the bound: w > 4=. We also find 
that there is no local modulus minimum for large axion 
field values: 



> 



V6 



no modulus minimum . 



(13) 



That is to say: for large axion field values, the corre- 
sponding modulus will become destabilised^ . The disap- 
pearance of the minimum at large axion field values is 
precisely the appearance of the "valley walls" in the po- 
tential, as described in [49] . 

When the axion has a periodic potential, canonically 
of the form U{0) — 1 — cos(6'), then the axion has a 
maximum field value at = tt. Such a periodic field can 
spoil the local modulus minimum when: 



D 2 



(14) 



Eqs. (13 14 ) show that if the natural scales in the axion 
potential (cither or /i^) arising from non-perturbative 
physics, are of the same order or slightly larger than the 
natural scales in the modulus potential (in this case £*), 



which are also non-perturbative, then destabilisation can 
occur even for small field values. Whether or not this 
mild hierarchy of scales occurs in actual models of moduli 
stabilisation is not the subject of this work, but we see 
no a priori reason why it should not be possible. 

If the axion initial conditions are such that field values 
are large then the corresponding modulus has no poten- 
tial minimum in the early universe. If the axion under- 
goes monodromy [551 - 155] the shift symmetry is broken 
and large field values are natural. For smaller axion field 
values this condition can still be satisfied for sufficiently 
large /3^, or if the bound of Eq. (14) is satisfied. 



In the rest of this section we will be concerned with 
situations where a local minimum for x does exist. In this 
case where there is a local modulus minimum, there is 
still interesting physics caused by the axion background. 
If the modulus begins life at its local minimum in the 
frozen axion background then the fractional change in 
the modulus field during axion evolution from (j) = (pi to 
(i = is: 



Ax I A/a I 



In 1 - 



(X>(0) fa,^ 



(15) 



This ratio blows up when <f>i = cf), where the modulus is 
destabilised and the local minimum is at ^ = 0. How- 
ever, it remains 0(1) for </> < (j)^ 1 — As we will 
see, Lo does not appear in our dynamical system analysis 
and so can be picked arbitrarily large (corresponding to 
stabilising the modulus at larger and larger values) and 
these results can be made insensitive to the approxima- 
tion that fa is fixed. For the consistency of our assump- 
tion that PQ symmetry is broken before inflation, we 
must have fa larger than the inflationary energy scale, 
and it must remain large enough that the symmetry is 
never restored. Two comments are in order here. Firstly, 
the modulus will never roll out to truly infinite values 
because of Hubble friction. Secondly, however, we may 
in general expect trajectories that go from a destabilised 
region on the plateau of the potential into the global min- 
imum to require some fine tuning in order not to break 
our assumption of fixed fa- We will comment more on 
this later. 

The ratio of modulus mass at the start and end of this 
trajectory is: 



™x((x)(0»)) 
™x((x)(0)) 



(16) 



This ratio can become small as the bound of Eq. ( 13 ) 
becomes saturated. When this bound is saturated, or 
nearly saturated, and if other corrections to the modulus 
mass are small^, then an effective field theory obtained 



® The same is true when the variation of fa with x is taken into 
account, but the condition must be found numerically. 



For example, suitable decoupling occurs in the scenario of [9]. 
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by integrating out the modulus based on its mass at the 
global minimum may fail. The potential is anharmonic, 
and so if the modulus is displaced far from its local min- 
imum the mass will not be given by this formula. 



The ratio of Eq. ( 16 ) occurs also for the axion mass 



along this trajectory, implying that in such a situation 
the fields cannot change their relative masses during the 
course of their evolution. If the modulus is heavier at 
the global minimum, it will also be heavier in any local 
minimum. This allows for consistency of the assumption 
above: if the modulus minimises first it should be a good 
approximation in this case to consider the trajectory as 
being We will call this trajectory the adiabatic 

trajectory^ i.e. the one that the modulus follows if it is 
always heavy enough to return to equilibrium suHiciently 
quickly as the axion rolls. 

We can obtain an effective potential for the axion that 
approximates the full effective field theory description 
by substituting the adiabatic trajectory for the modulus 
back into the potential of Eq. Q : 



2uj 



D 



6 



2uj 

^ PA 
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(17) 



This potential differs from the harmonic potential by be- 
coming flat as — > (/). Beyond (j) = <j> it should not be 
used. 

During this evolution the fields follow a curved tra- 
jectory in field space, just like in multi-field infiation, 
with the axion and modulus both always moving to 
smaller values and becoming heavier. Therefore the nor- 
mal course of cosmic evolution will not endanger late 
time stability. However, as the bound of Eq. (13) be- 



comes saturated we should see that axions and moduli 
undergo significant evolution in their masses while mov- 
ing towards the global minimum. If the axion is light in 
the current epoch then this evolution will still be occur- 
ring. 

Again drawing the analogy to multi-field inflation, even 
if the evolution in the modulus direction is slight, a tight 
turn in the fleld space trajectory may lead to observable 
features in the axion power spectrum j86| I87j . Such tight 
turns do not appear possible in the potential we study, 
and we also will not be considering the effect of inhomo- 
geneous perturbations. 

We finally note here that if the amplitude of axion 
oscillations were allowed to grow, such as in the scenario 
explored in |25j where axion oscillations are amplified 
by the Penrose process near a black hole, then this may 
also lead to novel effects on the modulus sector and the 
vacuum energy. 



Dynamics in Phase Space and the Equation of 
State 



The value of the vacuum energy (which should not be 
confused with the scale of the potential in Section ([lT|)), 
A = V{0,{x)), which includes the pA contribution, is a 
free parameter in our model. Its sign controls the topol- 
ogy of phase space [88, . We show this effect in our model 
in Fig. ([2]), where we plot trajectories in {(j>, (j), H} phase 
space obtained by numerically solving the equations of 
motion (Eqs. ([6|). With A > the trajectories are con- 
fined to the expanding branch, H > (or if Hi < 0, the 
contracting branch), and the phase space is disconnected. 
With A < it is possible for the total energy density to 
go to zero, and so — > connecting the expanding and 
contracting branches and making phase space connected. 

When H = the evolution of the scale factor turns 
over, such that with A < the universe expands and then 
contracts to a big crunch despite, in these cases, having 
zero curvature (see e.g. (HH] and references therein). We 
plot the evolution of the scale factor for A < in Fig. (|3|. 

The trajectories of Fig. ^ all begin at the local mod- 
ulus minimum, {x){4>i), with stationary fields in a fiuid 
dominated universe, and (pi takes various values between 

and (j). We begin at time t — Q^a — 1 during matter 
domination and look at the evolution towards the dark 
energy universe of today. We have shown trajectories 
given by evolution in the full potential (dashed lines), and 
in the effective potential (solid lines). Those trajectories 
with large initial axion field values in the full potential 
depart from the evolution in the effective potential. This 
is because at large axion values the modulus is light and 
the adiabatic assumption is no longer good enough. The 
scale factor evolution for A < is shown in Fig. ([s]) and 
a clear difference is visible between evolution in the two 
potentials, with the maximum size of the universe being 
larger when the full potential is used. 

We investigate the accuracy of the adiabatic approxi- 
mation for A > in Fig. Q, where we give Ax/x(0 = 

1 — {x){4>)lx in percent. We see that at early times the 
trajectories with large initial (f) depart by as much as 
25% from the adiabatic trajectory. All trajectories un- 
dergo damped oscillations about the adiabatic trajectory 
as the modulus mass increases over time. 

We conclude this section by commenting on the ef- 
fect of the combined axion-modulus oscillations about 
the minimum on the dark energy equation of state. With 
A > and at least one light field (conservatively, an axion 
of the axiverse) , the fields at or close to their initial values 
can come to dominate the energy density as dark energy, 
yet still be evolving towards, and oscillating about, the 
true vacuum at late times. We show the effect of this on 
the equation of state W0xA — {p^+p^+piC) / {P<t,+Px— Pa) 
in Fig. ([5]), for the trajectories of Fig. ([2| (left panel), 
again comparing the cases of the effective potential and 
the full potential. 

In this evolution, the dark energy axion-modulus fluid 
is already dominating the energy density at early times. 



FIG. 2: Phase space topology in {(j), (j), H}. Left panel: A > 0, with evolution to a minimum H. Right Panel: A < 0, trajectories 
spiral through H — and the universe collapses. Dashed lines are for evolution in the full potential, and solid lines in the 
adiabatic potential, where the modulus remains always in its local minimum. 
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t 

FIG. 3: The evolution of the scale factor for A < 0, showing 
a turn over and collapse of the universe despite there being 
zero curvature. Again, dotted lines are for evolution in the 
full potential, while solid lines are in the adiabatic effective 
potential. Bottom to top corresponds to increasing — > (j). 
At large (f) the universe reaches a larger size when the full 
potential is used. 




FIG. 4; Comparing modulus evolution for motion in the full 
potential versus motion in the adiabatic effective potential 
with A > 0. The modulus always begins in the local mini- 
mum. Lines from bottom (light blue) to top (dark blue) repre- 
sent increasingly large axion initial field values. The modulus 
departs by up to 25% from the adiabatic trajectory when the 
initial value of the axion field is large. 



t ~ 100, with w « — 1. However, w is rising as the fields 
move towards the minimum leading to large departures 
from w = — 1. For large axion initial values in the full 
potential this motion is delayed and w remains flatter for 
longer. In all cases, as oscillations begin w rises so much 
as to halt accelerated expansion altogether {w > —1/3). 
The positive value of A in the true vacuum means that 
at late times w will relax to exactly —1 and that w is 
bounded to — 1 < w < 1. This bound does not hold for a 
negative potential [SI], indeed > 1 was observed for 



this potential in 



IV. DYNAMICAL SYSTEMS ANALYSIS 

To delineate the regions of parameter space that may 
give rise to acceptable and interesting cosmological phe- 
nomenology we perform a dynamical systems analysis by 
transforming the coupled axion-modulus system into au- 
tonomous form. The axion and modulus fields and the 
perfect baryotropic fluid evolve according to Eqs. (|6|, 
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FIG. 5: The Dark Energy (DE) equation of state for A > 0. 
Dotted lines are for motion in the fuU potential, solid lines 
for motion in the adiabatic effective potential. Left to right 
lines have increasing — > 0. DE has come to dominate the 
energy density around t — 100, but later both the DE and 
total equation of state go back above w — —1/3 (horizontal 
line), temporarily halting accelerated expansion. In the most 
extreme case of large 0i in the full potential this can happen 
even when the equation of state today is very flat and close 
to w = —1. The late time expansion is asymptotically dS as 
the fields relax into the minimum. 



subject to the Friedmann constraint Eq. ([t]). 



with 
H' 



y = 



t' = ^ 



H' 
ll 

IT 
ll 

ir_ 
li 

IT 

ir_ 

H' 



X — Prs , 
3 



y- 

3 



C\r^ - - 2z^] 



Cy 
Cy 



Pxr , 



H 



(19) 



-'fb{l-x^-s^-y'^-z^+r^-t^)-ix^-iy^ , (20) 



where (3 was defined in Eq. (12). 



Here, a prime denotes diiTerentiation with respect to 
the number of e-foldings N = In (a). The dimensionless 
density parameters fl^ = pi/SH"^ of the cosmic compo- 
nents i can be expressed as 

nA = t\ (21) 

t^) . (22) 



— x'^ + s'^ , 11^ = y~ + z~ — r~ 
and furthermore, fiatness imposes 



rJb = 1 - (a:^ + + + -r^ 



A. Autonomous System 



Following [3U] we make the change of variables: 




t = 




(18) 



The evolution Eqs. ^ can then be transformed into au- 
tonomous form X' = f(X), where X is the column vec- 
tor of compact variables and f(X) is the corresponding 
column vector constituting the autonomous system equa- 



At this point some comments on the system (19) are 
in order. Notice that due to the negative contribution 
from Vd in the modulus potential, trajectories are not 
confined to the unit hypersphere in the full phase space. 
The set of phase space variables {x, y, z, r, s, t} is of one 
dimension more than the actual {(j), (j), x, x, H} degrees 
of freedom. This is because there is a relation that ex- 
ists between the phase space variables, which provides an 
additional constraint and defines a surface on which the 
motion takes place, just like the Friedmann constraint 
gives the topology of phase space in [88] . The constraint 
is: 



t^z^ 



BpA 



(23) 



which is a simple consequence of the definitions of the 
variables and the form of the potential. Trajectories are 
confined to live on this plane, defined by the choice of ini- 
tial conditions. There is a simple interpretation of this 
that will help us visualise phase space: Choosing ^ corre- 
sponds to a choice of sign for the vacuum energy, includ- 
ing the bare cosmological constant. Combining Eqs. ^ 
and ( 11 ) and ip — to get the vacuum energy {V) ~ A: 



B V A 



(24) 



10 









Zc 




Sc 


tc 


Existence 


A 




















all p,C,7b 


B 

















±1 


all C 7h 


C 




±y 














-1 <y< 1 


D 


















\/f 


E 





V 8 C 













7b < 2 


F 





C 








V 6 





C < \/6 


G 





pi lb 

y 2 c 













7b < 2 


I 





C 
V6 













c>V& 


M 








±z 


±^/2z 





±Vz2 + 1 


all /3,C,7b 



TABLE 11: The fixed points of the system (191 and the conditions for their existence. Rather than having an isolated fixed 



point, M is formed of a continuous line of fixed points, called a critical line. This critical line intersects the C, plane at a unique 



point z = zm given by Eq. 1 30 1 



The Minkowski vacuum is given by C = whilst C, < \ 
is an AdS vacuum, and ^ > | is a dS vacuum (we show 
the planes for three values of ^ in Fig. (|6|, which we will 
discuss in more detail later) , but we note that C, does not 
set the scale of the vacuum energy. This would be fixed 
observationally if we were to want the vacuum to give 
us the correct iJp, but this is not necessary: {V) can be 
much less than this, but not greater. 

Furthermore, this geometrical picture of surfaces in the 
{z, r, t} subspace can give another view on the bounded- 
ness of trajectories and topology of phase space discussed 
in Section (III) in Fig. The requirement of a flat uni- 



verse imposes the constraint: 



V 



2 2 

z — r 



<1. 



(25) 



Clearly, the location of the hypersurface defined by sat- 
uration of the bound moves as the variables evolve, how- 
ever we can picture its effect in the limit of heading to 
the vacuum: x = y = s = Q. Now there is an additional 
surface that intersects those of Fig. ([6]). If one were to 
plot it, one would sec that it intersects dS surfaces, with 
C > 1/4, making an arc below which trajectories are con- 
fined, unable to reach co-ordinate infinity. For AdS sur- 
faces, with C < 1/4, the surface funnels outwards, re- 
stricting trajectories to a region of their C surface, but 
not confining them to finite values. This is another man- 
ifestation of our choice of dynamical system variables: for 
a negative vacuum energy it is possible for iJ — >■ 0, where 
the co-ordinates diverge and trajectories on the potential 
become unconfined. 

The plane defined by C = j, which we will call 
the Minkowski plane, therefore divides the phase space 
into three: above the Minkowski plane, phase space 
is bounded, and the collapsing and expanding universe 
branches (z, r,s,t < and z,r,s,t > respectively) are 
separated; on the Minkowski plane the asymptotic fu- 
ture lies at infinity, where H = 0; below the Minkowski 
plane, phase space is connected and there are trajectories 
through infinity that join the expanding and contract- 




FIG. 6: Dynamical constraint surfaces in the {z, r, t} sub- 
space. The fiat meshed (blue), surface corresponds to the 
Z plane, which is the minimum of the potential. The fore- 
most curved (yellow) surface is the constraint for a dS vacuum 
with > 1/4. Next behind it is the Minkowski plane (green) 
with ^ — 1/4. Finally, the rearmost curved surface (red) is 
the constraint for an AdS vacuum with ( < 1/4. The bold 
(red) line on the Z surface corresponds to the critical line M. 
We see that M crosses only dS planes and asymptotes to the 
Minkowski plane at co-ordinate infinity, while the minimum 
surface crosses all ^ planes on a line. The non-crossing of fixed 
line M with surfaces of < 1/4 is another geometric illustra- 
tion of the instability of a negative potential to collapse. 



ing branches. When the transition is made from expan- 
sion to contraction and H changes sign, all six variables, 
x, y, z, r, s, t also change sign. 

Choosing to work with a phase space of one higher 
dimension is useful for our analysis of fixed points, since 
it allows us to see what happens when one or more of 
these variables can be approximated as vanishing relative 
to the others, for example z and r vanish as the modulus 
goes to large values, and t vanishes in the early universe. 
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TABLE III: Properties of the fixed points given in Tabic (|ll| for an expanding universe. For M, the contribution from pA has 
been included in Wcff. 



When H — our variables diverge and so the system ( 19 ) 



cannot be evolved through the transition between ex- 
panding {H'^) and contracting {H~) universes. It is ac- 
tually possible to construct a set of compact variables 
which remain finite at H — 0: 



V2Q ' 



VQ 



X 

V2Q ' 

Q ' 



/PA 



Q 



where 



Q = ^3H^ + Vd 



(27) 



These variables are similar to those defined in ^T] . Since 
Vd is positive definite, Q always remains well defined. 
Defining a new independent variable (') = one 
can transform the evolution Eqns. ^ and Q into au- 
tonomous form. This alternative autonomous system is 
given in Appendix ([b|. These compact variables, Xq, are 
related to our original compact variables X (Eqs. (18)) 
by 



Xqf 



(28) 



where r was defined in Eqs. (18). We have kept our 
variables finite at H = 0, at the expense of losing 
the intuitive description of the division of phase space 



provided by the vacuum constraint Eq. (23), since the 



Xq system has the minimally required dimensionality. 
For this reason we content ourselves with describing 
the axion-modulus system in terms of the variables of 
Eqs. (18), and do not study the transition at H = 
explicitly. We will occasionally make use of the Xq 
variables to numerically show the evolution of phase 
space trajectories. 



B. Fixed Points 

The fixed (critical) points X^ of the autonomous sys- 
tem ( 19 1 are extracted by satisfying X' = and are 



listed, along with their conditions for existence, in Ta- 
ble (|ll]). As mentioned earlier, the positive (negative) 
roots in the {z,r,s,t} subspace correspond to expand- 
ing (contracting) universes. In total there are thirteen 
fixed points, four of which are imaginary and so are not 
physical and are not listed in Table (|ll|. The energy 
densities fii, the effective scalar field equation of state, 
(26) ""^cff = (^V+^x)/(P0+Px) conditions for acceleration 
and stability of these fixed points are given in are given in 
Table (III). A fixed point corresponds to an accelerating 
solution if 



(1 - x2 - yl - si - zl + rl - tl)-fk + 2x1 + "^vl < 3 ■ (29) 

The stability of the fixed points may be determined 
by expanding about them, setting X = Xc + (5X, with 
(5X the perturbations of the compact variables defined 
by Eqs. (18) considered as a column vector. To first or- 
der, the perturbations satisfy 6X.' = W • JX, where the 
matrix W contains the coefficients of the perturbation 
equations. The stability of the fixed points thus depends 
upon the nature of the eigenvalues of the matrix W. The 
full stability analysis is somewhat cumbersome and may 
be found in Appendix ( [A| . Here, we give a general sum- 
mary of the fixed points and their stability, focussing 
on the intuitive physics that the dynamical systems ap- 
proach provides. 

Of the nine fixed points listed in Table (|ll]), there are 
two trivial solutions: Fixed point A corresponds to the 
fiuid dominated point where the kinetic and potential 
components of the axion and modulus fields are negli- 
gible, whilst fixed point B represents the p\ dominated 
solution. Point A is unstable in both expanding and 
contracting universes. Recall the ultimate fate of the uni- 
verse is determined by the value of ^. In the presence of 
a dS vacuum {C, > 1/4) the stability analysis reveals that 
fixed point B is associated with three zero eigenvalues 



12 



in the {z, r, s} subspace, whilst the remaining directions 
are stable. We say that this is a marginally stable so- 
lution in the sense that there is no instability growing 
exponentially, although it could be unstable to higher 
orders in the perturbation. To obtain the strict stability 
of this solution we would have to go beyond linear order 
in perturbation theory, which we do not pursue as nu- 
merical integration of the autonomous system confirms 
that this point is ultimately unstable: the asymptotic fu- 
ture in the presence of a dS vacuum is the stable fixed 
point M, the global axion-modulus potential minimum, 
which has a larger basin of attraction. The existence of 
point B demonstrates the ability of a bare cosmological 
constant to overdamp modulus motion for the modulus 
beginning life high up on the plateau of its potential, 
shielding us from the true vacuum and seeing only the 
larger pA- In Fig- ^ we show this temporary 'trapping' 
in fixed point B by plotting trajectories in the {zq^Iq} 
subspace. Whilst such a trapping may last for hundreds 
or even thousands of e-foldings, the modulus will even- 
tually begin to roll when its mass overcomes the Hubble 
damping and will relax into its minimum. 

Rather than having an isolated fixed point, point M is 
formed of a continuous line of fixed points, known as an 
equilibrium manifold, which we call a critical line. The 
emergence of this critical line is due to the fact that we 
are working in one dimension more than is required. In 
one dimension less, the line would degenerate to a unique 
point, which is given by the intersection of M with the 
(■-plane. 



ZM = ± 



1 



(30) 



which clearly only exists for ( > -g- That is to say, the 
absolute potential minimum defined by the stable fixed 
point M is only a fixed point in the presence of a dS vac- 
uum and corresponds to the asymptotic future. If C < j 
the global minimum is not a fixed point, and the asymp- 
totic future is cosmic doomsday in a Big Crunch as will 
be discussed below. The global minimum is expressed in 
terms of the autonomous system variables by substitut- 
ing z = zm in M. We finally note that the line M can 
equally be derived as the intersection of the minimum 
surface with the saturation of the flatness constraint in 
the potential dominated regime, z^ — r"^ + = 1. 

As well as confirming that M is stable for C > 5, the 
stability analysis reveals two bifurcation points that lie 
along it: 



3 
8C2 



and 



9 



(31) 



These bifurcation points are obtained by setting the 
quantity under the square root in the eigenvalues Ai_2 and 
A4_5 of Eq. (All ) to zero and solving for z. For M to be a 
stable node in the {y, z, r} subspace, z < y/3/8C^, other- 
wise it is a stable spiral, whilst for point M to be a stable 
node in the {x, s} subspace, z < -y/O/S/?^, otherwise it is 





a stable spiral. The vacuum surface C, determined by B, 
Pa and D, dictates which side of the bifurcation points 
Zm lies. We have the conditions 



(32) 



which are derived by setting z = zm in Eq. (31 1. The 
stability analysis has elegantly revealed the conditions 
for late-time oscillations of the axion and modulus fields: 
violation of the /3 condition corresponds to axion oscilla- 
tions, whilst violation of the C condition corresponds to 
modulus oscillations. Fig. ([7| shows an example of late- 
time modulus oscillations as the trajectories spiral into 
the point zm- The time scale of the axion oscillations 
are determined by /3, C and also If the vacuum is 



chosen to be Minkowski, the conditions (32) simply be- 



come /3 < and C < 0, which are never satisfied for the 
parameter values considered in this work and so axion 
and modulus oscillations are inevitable. 

There is another plane which is of interest in the 
{z, r, t} subspace, which defines the minimum of the po- 
tential itself at (x). At (x), = V^z, which defines the 
plane, and the critical line M lives here. This plane, 
which we will call Z, crosses the ( plane on a line. Tra- 
jectories along this line are those living in the minimum 
and leading to the asymptotic future, either at the cross- 
ing point of M in a dS vacuum, or ultimately leading 
to collapse in an AdS vacuum. Trajectories crossing this 
line are modulus passages through, or oscillations about, 
the minimum. Trajectories in the full 6-d space, how- 
ever, never cross each other: these are oscillations and 
static passages along the crossing and are separated 
in the y-direction. These surfaces are shown in the ex- 
panding octant of the {z,r,t} plot in Fig. (|6|, where we 
see the crossing of Z along a line in the dS, Minkowski 
and AdS example C. planes, M lying in the Z plane, and 
crossing the dS plane at a point. 

Fixed point C is the second critical line of the system, 
corresponding to an axion-modulus kinetic dominated 
(stiff fluid) solution. This critical line is the unit circle 
x^ + y^ = 1 and is a symmetry of the autonomous system 
with Zc = rc = Sc = tc = 0. This is the usual enhance- 
ment of symmetry for massless scalar fields. In an ex- 
panding universe (where stability is ensured by negative 
eigenvalues of W), this point is always unstable. For a 
collapsing universe, a fixed point is stable if the eigenval- 
ues of W are positive. This is because the 'time' variable 

= In (a) of the autonomous system becomes a decreas- 
ing function of time. Hence, critical line C is stable in 
a collapsing universe and corresponds to the asymptotic 
future of any model with an AdS vacuum. This is con- 
sistent with the pre-big bang cosmology j2l [13] late time 
attractor solutions. The particular fixed point along C 
that the system will finally evolve to will depend upon 
the initial conditions of the system. Similarly to the 
phase space dynamics in the presence of a dS vacuum, 
the only possibility to save us from this Big Crunch cos- 
mic doomsday is a temporary trapping in fixed point B. 
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FIG. 7: The evolution of phase space trajectories in the presence of a dS vacuum ((" = 0.276) obtained by integrating Eqs. ( |B1[ |. 
We set C = 1, 7b = 1 and /3 = 1.87. The compact variables zq and Iq are related to z and t through Eq. (281. Left panel: The 



temporary trapping of the zq and Iq trajectories in fixed point B b efore the modulus begins to roll, finding its minimum at 



fixed point M. Saturation of the bound C < | (i^ — |) , Eqn. (321, results in late-time modulus oscillations, which are seen 
in the figure as trajectories spiralling into M. Right panel: tq as a function of TV, the number of e-foldings. 
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FIG. 8: The evolution of phase space trajectories in the presence of a dS vacuum obtained by integrating Eqs. ( 19 1. Left Panel: 
Trajectories repelled from the unstable fixed point F and heading for the global potential minimum at M. We set ^ = 1.0, 
C = 0.5, 7b = 4/3 and /3 = 1.732. Right Panel: Trajectories spiralling into the unstable fixed point G, before the modulus 
finds the global potential minimum at M. We set C, — 1.0, C = 10.0, 7b = 4/3 and P = 10.0. 



This situation was seen in the examples of |49j whenever 
Pa domination sets in before collapse and is achieved for 
large initial modulus values. This period of dS inflation 
would only be temporary however and the modulus will 
quickly relax into its AdS vacuum signalling cosmolog- 
ical collapse and leading to eventual decompactification 
as described in This is demonstrated in the exam- 
ple plot of Fig. ([9): the growing kinetic energy of the 
modulus as a — ?> drives it to large values. 

All other fixed points are unstable in the presence of 
an AdS or dS vacuum. Points D and E correspond to dy- 
namical modulus stabilisation at small modulus values, 
while I corresponds to dynamical stabilisation at large 
field values. This is of course only a meta-stability, since 
these fixed points are unstable. Point E is also a scal- 
ing solution, on which the axion energy density vanishes 
(p^ = 0) and the modulus energy density scales with the 
dominant background fluid: 



Px 



9Hf 

4C2 \ a, 



-37b 



3 7b (7b 
4 



1) 



(J2 



(33) 



Hence, the modulus tracks the dominant background 
fluid and Px/ Pb remains constant. 



Fixed point F represents a solution dominated by the 
modulus kinetic energy and the potential energy of the 
axion. On this solution, their relative energy densities 
and effective equation of state remains fixed: 



Px 



6 



1. 



1 



1. 



(34) 



The repulsive nature of fixed point F is illustrated in the 
left panel of Fig. Q. The only fixed point which admits 
a non-vanishing background fluid density with a sizeable 
contribution from both axion and modulus is G. Here, 
the modulus energy density is dominated by its kinetic 
contribution, whist the axion remains frozen, its motion 
suppressed by Hubble friction. Both the axion and mod- 
ulus track the evolution of the dominant background fluid 



^Hf .fa 



-37b 



Px 



2C2 v a. 



-37b 



(35) 



whilst giving a background density J7b = 1 — "i^hjC^- 
This dynamical attractor is precisely the axion and mod- 
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FIG. 9: Evolution of the modulus field in a collapsing uni- 
verse. As a — ^ at t ~ 10~^ the kinetic energy grows and 
dominates in fixed point C and the field value diverges. Since 
this would happen to all scalar fields and hence all moduli, 
this signals decompactification near a crunch. In this example 
plot the units and parameter values are all arbitrary. 



ulus tracking behaviour that was described in |49]. It is 
in the combined equation of state 



WcS 



37b (7b 



1) 



(J2 



(36) 



(rather than the individual equations of state) that we 
see tracking as Wgs tries to follow the equation of state of 
the dominant component. Tracking is finally destroyed 
as axion oscillations begin, which is the cosmic trigger 
event that restabilises the modulus. We show evolution 
into this fixed point in the right panel of Fig. (Isl). 



Additional Comments on Fixed Point G 

As a particularly interesting fixed point, we choose to 
discuss some additional phenomenology relating to fixed 
point G. Firstly we discuss accessibility of the fixed 
point. Even in the unbounded co-ordinates of an AdS 
minimum, approximate trapping in G requires variables 
other than {y, s} to be approximately zero, and so flat- 
ness bounds us with + < 1, which defines a circle. 
This in turn imposes a constraint on C as a function of 
7b for G to be within this region: 



C > 



physically accessible G. 



(37) 



An interesting phenomenon when entering fixed points 
in a multi-field model is the possibility of multiple periods 
of accelerated expansion [93]. When projected down to 
the {y,s} subspace the condition for w < —1/3 in an 

)• 
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expanding universe becomes s > y ( ^^'' ){y 

A temporary trapping in point B, where there is a 
larger value of cosmological constant than in the true 
vacuum, could lead to a single period of accelerated ex- 
pansion during an otherwise radiation dominated era. In 
such a case, the axion and modulus fields would pick up 
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FIG. 10: Constraints relevant to trajectories approaching the 
quasi-stable fixed point G with a background matter fiuid, 
7i, = 1. The blue shaded region under the semi-circle is the 
region allowed by the Friedmann constraint. The hyperbola 
bounds accelerated expansion, with the red shaded region 
having w > —1/3. The two points represent fixed point G 
for s > 0. For C > \/3 these lie inside the allowed region. 
We conclude that it is possible for trajectories approaching 
G to cross the w = —1/3 divide if they spiral as they do so, 
possibly leading to multiple epochs of accelerated expansion. 



large additional isocurvature fluctuations from this brief 
period of inflation. This period would end when the fields 
move towards their vacuum, and as such the global min- 
imum would have to be dS. 

The situation for G in this regard is more interest- 
ing. In a matter background, — 1, this is pictured in 
Fig. 10 Here we see that it is possible to have both a 
fiat universe and w < —1/3 in a {y, s} dominated phase. 



Trajectories in the {y, s} subspace will spiral into G if 
the eigenvalues of the stability matrix W (that point in 
the {y, s} directions) have an imaginary part. The sta- 
bility analysis (see Appendix (|A|, Eq. ( |A8| )) reveals that 
this is the case if: 



- SC^ + 24 > . 



(38) 



For the parameter space of interest, C > 0, this bound 
is satisfied for C < 2 and C > 3 + ^/2i « 7.6. For any 
C between these two values, the eigenvalues are real and 
trajectories will not spiral into G but move in straight 
lines, and so cannot cross w < —1/3. When Eq. (38) 



is satisfied however, the trajectories in {y, s} can spiral 
toward G, having the possibility of crossing the w < 
— 1/3 bound, perhaps multiple times. So, trajectories 
approaching this fixed point can lead to multiple periods of 
accelerated expansion during a matter dominated epoch. 
This phenomenon is extremely tightly constrained: such 
an epoch of acceleration must be less than 0.05 e-folds 
long [M]. 
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C. Scanning Parameter Space 



Since the system of Eqs. ( 19 1 are first order and au- 



tonomous, they are very quick to integrate numerically. 
We exploit this nice property by performing a 'scan' of 
the model parameter space around regions of interest, 
selecting particular scenarios to investigate more system- 
atically. We use our scans to further our qualitative un- 
derstanding of the phenomenology of the model and to 
locate and single out specific novel features. 

The autonomous system has eight different param- 
eters which determine the subsequent motion of any 
given trajectory in phase space: six initial conditions, 
{xi,yi, Zi,ri, Si,ti} and two parameters, {C, /3}. To en- 
sure that this rather large parameter space is sampled in 
a uniform and efficient way, we use the method outlined 
in Appendix ([C|). 

We briefly describe this process for initial conditions 
chosen to be close to fixed point A, i.e. beginning in the 
fluid dominated phase with a non-vanishing background 
fluid density, ilb (initial). It is trivially generalised to the 
case of any other fixed point. For A, with only some loss 
of generality, we make the simplifying assumption that 
the axion and modulus fields begin frozen, Xi ^ = 0. 



Then, using the Friedmann constraint, Eq. (22), and the 
vacuum constraint, Eq. (23 1, we have 



(39) 



initially. Here, p = 1 — (initial). The initial conditions 
are constrained to lie on this three-dimensional manifold, 
which we will call Ai. Scanning the initial conditions of 
the system then reduces to varying two initial conditions 
evenly over M with the third constrained by the equation 
for M. We choose to vary Zi and r^, whilst still being free 
to independently vary {C, /3}. Since M has non-constant 
curvature, it is not trivial to sample it in a uniform way 
and so we use a statistical sampling method which is 
presented in Appendix ([C]). Choosing a value of p ~ 0.01 
(fib (initial) « 0.99) is our definition of "near" to fixed 
point A. 

Before we present the results of our numerical scans, 
it will further add to our intuition to briefly discuss the 
change in topology of the initial condition manifold either 
side of C = 1/4. If the vacuum is dS, C > |, the surface 
area of the manifold above some value s — Smin is finite. 
For C < J, (AdS and Minkowski vacua) the manifold is 
not bounded above Smin a-nd its surface area is infinite. 
Suppressing the subscript i for brevity, this change in 
topology is best illustrated by solving Eq. (39) at s = 
Smin for z. This generates two physically relevant roots 
as a function of r which describe the curves where A4 
intersects the s ~ Smin plane. These two curves meet at 
a point, r-max, 



One can also obtain equations for the two A^-Smin inter- 
section curves as a function of z. These two curves meet 
at 



UCiP- 



4C-1 



(41) 



Eqs. (40) and (41) illustrate the change in topology of 
M: for C = 1/4, r^ax , ^max oo, whilst for C < 1/4, 
''max , -Zmax G C. In both cases, the manifold never 
intersects the s = s^in plane. Only for ^ > 1/4 is the 
surface area of the manifold bounded above Smin- This 
is another clear example of the genuine change in the 
geometry of phase space when the vacuum is chosen to 
be either dS or AdS. 



We now begin to discuss the results and findings of 
our numerical analysis. Our ability to perform scans of 
this kind has many possible applications for investigat- 
ing the cosmological phenomenology of our model. Here 
we choose to simply show some examples that illustrate 
the capabilities of our technique. All models we present 
have a dust background fluid, 7b = 1 and a dS vacuum. 
We run two large simulations: FP-A and FP-G. For sim- 
ulation FP-A, we evolve 562500 models, each chosen to 
begin in fixed point A with ilb(initial) 0.99. We scan 
{C, f3} evenly in logarithmic space on a 25 x 25 grid, and 
at each point we use our initial condition algorithm to 
evenly sample the space of {zi,ri,Si} over M. For FP-G, 
we evolve 506100 models, each chosen to begin near to 
fixed point G. Unlike simulation FP-A, we are not free to 
independently vary the initial conditions and C since in 

fixed point G, Ob = 1 — 3/C^. Furthermore, j/c — \J~^/h 
and so yi ^ 0. Hence, every time C and yi are changed, 
the shape of the initial condition manifold also changes. 
Ther efor e we absorb yi and Ob into the parameter p of 
Eq. ([39]): 
cally and 

initial condition algorithm to evenly sample the space of 
{zi,ri,Si} over M for each point, {C, j/i}, where Ai has 
a different shape. /? is varied 15 times on a logarithmic 
scale. 

In both simulations, individual models are terminated 
under two conditions: either they have settled into fixed 
point M for more than 5 e-folds, or, they have ran for a 
total of more than 500 e-folds. The results are presented 
so that at each point in {C, /?} space, the average over all 
trajectories on Ai is taken, or alternatively for each point 
in {zi, Ti} space we could average over parameters {C, /?}, 
i.e. repeated points in any plane have their contour value 
averaged. 



p = ^ — yf- We then vary C logarithmi- 
yi linearly across a 10 x 15 grid and use our 



The End of Fluid Domination 



' (p-^Ln)(l+2^/C) 

4C-1 



(40) 



In Fig. (11) we plot in the {C, /3} plane, for models 



from the FP-A simulation, the number of e-folds, iVofd, 
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FIG. 11: Contour plot of the number of e-folds, A'ofd, when 
Q,b < 0.5 for a scan over parameter space of all models be- 
ginning near fixed point A with Qb ~ 0.99 (the FP-A simu- 
lation). Each point in {C, /?} represents an average over the 
initial condition manifold A4. 



When considering models from FP-G we imposed a cut 
for all C < VS, where G is unphysical (the initial con- 
ditions correspond to negative Jib). The first difference 
observed from fixed point A was vertical cut giving very 
low A'ofd at small C. These models had initially very 
small fit: the minimum of A'ofd output from the code 
is A'ofd = 0.005 which is our numerical step size in A', 
i.e. the models began out of fluid domination. Since C 
sets the initial fit for models beginning in G, the gen- 
eral trend of increasing A'ofd with C continued and was 
dominant, until at large C and j3 it gave way to the ef- 
fects described above in the case of A. Axion oscillations 
lead to a decrease in (p to cj) < (j) and therefore spoiled G 
after some short time, decreasing the overall scalar field 
density. In the case of an AdS negative potential min- 
imum, these regions where G is spoiled would be those 
that eventually collapse. 



2. Multiple Periods of Accelerated Expansion 



from the beginning of the evolution until the end of fluid 
domination when fib < 0.5. 

Before analysing the figure in more detail, it is first 
worth thinking about what we should expect from such 
a number. The maximum across all models occurred for 
A^efd ~ 5.3 while the minimum occurred for A'ofd ~ 1-5. 
It is simple to show from the Friedmann equation that a 
ACDM cosmology beginning with Jim (initial) = 0.9 will 
reach — flA — 0.5 after A'ofd ~ 1-5 e-folds. This 
is because, in the absence of energy input, a cosmologi- 
cal constant maximally decreases (1 — flm)'- this should 
be the limiting case of our model when the fields are 
frozen, which indeed it is. On the other hand, a model 
beginning at matter-radiation equality at a^q ~ 10"'^ has 
(1 — Jim (initial)) ~ 0(10^^), and depending on o has 
6 < A^cfd < 7. 

In our model A'ofd can be increased and approach this 
limiting case in three ways. The fields can oscillate be- 
fore they overtake the fiuid density, they will then scale 
like matter, always remaining sub-dominant and the end 
of fiuid domination will be caused by the cosmological 
constant. Secondly, they could enter a scaling solution, 
where they also remain a fixed sub-dominant fraction of 
the energy density. Thirdly, they could roll to the min- 
imum of the potential, reducing the vev due to the neg- 
ative energy term in the modulus only part of the po- 
tential. We see that our maximum of A'ofd approaches 
the limiting case, being slightly below it as some time is 
taken for these dynamics to occur. 



In Fig. (11) there is a clear correlation of A'ofd with 



the parameters. Smaller average A'ofd occurs for low (3, 
where the axion mass is small preventing oscillations, and 
low C where the scalar field energy density in scaling 
solutions is large. Larger average A'ofd occurs for large /3 
and C where oscillations can occur earlier and the energy 
density in scaling solutions is smaller. 



Motivated by the fact that spiralling trajectories in 
phase space may generate multiple periods of accelerated 
expansion, we scan the model parameter space for this 
feature. We compute the number of periods of acceler- 
ated expansion, A/'ao, by counting the number of times 
Wtotai = '^i'Wifli < — I along the model trajectory. 
Here, i labels the axion and modulus fields and the dust 
fiuid. If one, or both scalar fields are oscillating about 
their minima, the averaged equations of state, and 
are used in the calculation of wtotai- The average taken is 
a moving average and is re-calculated every 0.005 e-folds 
as the trajectory advances in time. We define the onset 
of coherent oscillations as the time when the field veloc- 
ity {(f) or x) changes sign for the third time. This ensures 
that we do not average any heavily or critically damped 
oscillations. This definition is somewhat arbitrary and so 
we should expect that A/'ao may be sensitive to the defi- 
nition of the averaging process. Furthermore, for regions 
of parameter space where the fields are highly oscillatory 
(large C and /3), sampling the trajectory every 0.005 e- 
folds may not be frequent enough to accurately average 
a single oscillation. We also note that this sampling rate 
is one tenth of the length of a period of accelerated ex- 
pansion allowed by observation. Hence, computing A/'ao 
by taking the moving average of and is not always 
the observationally relevant procedure. 

With these limitations acknowledged, we consistently 
apply our definition of A/'ao to every single model in our 
simulations. From the FP-A simulation we found that of 
our 562500 models, 17668 had Wao > 1. Of these, 676 
models were terminated after 500 e-folds for not reaching 
fixed point M, so that the multiple A/'ao can be said to 
have definitely occurred near to a fixed point or the local 
minimum. We also found 54 models with the largest 
A/'ao = 8. Of the remaining models with A/'ao — 1, 176599 
were terminated for not reaching M after 500 e-folds, 
and were thus still on the potential plateau trapped in 
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B. A/lc may increase in future for these models, but the 
time scale is immense: situating them today, 500 e-folds 
gives At = AN/Ho 10^^"^^ years^. We stress that we 
are not proposing any measure or figure of merit for fine 
tuning in this model, and as such the specific number of 
models pertaining to each case does not have any (clear) 
meaning. 

When considering the the distribution of A/'ac against 
{C, /?} we took all models with A/'ac > 1 and averaged 
over A4 as described above. We found some large regions 
of parameter space with A/'ac — 1 over all of Ai . We also 
saw that there was a high density of large A/'ac at larger 
values of C and intermediate values of /?, with one clear 
peak. We show these locations schematically in Fig. ( 12 ) 



Our results also showed an interesting correlation be- 
tween three dependent (output) variables where it was 
noticed that trajectories with large A/'ac occurred in those 
cosmologies that at the exit from fluid domination (en- 
tering the current epoch) contained only small values of 



ilrh and irJx 



(it is consistent in this model to have $7^ < 



since it does not contain p\: the total energy density re- 
mains always positive). This, combined with the larger 
values of C in these regions, as we will discuss below, sug- 
gests that these models were likely in or near to G (or B) 
at this time (again, see the schematic Fig. (12|). Small 
values of \rt^\ and fi^ for light axions are those allowed 
by current data (we discuss some bounds in Section ) , but 
is also potentially detectable with next generation exper- 
iments [ini IM] • Our scan suggests that such a cosmology 
could reasonably expect to undergo multiple periods of 
accelerated expansion in the future, and may have in its 
past. We re-state the bound from above: |S4J showed 
that an intermediate epoch of accelerated expansion in 
the matter era must have lasted less than 0.05 e-folds. 

Finally, our results showed that that almost all models 
with A/'ac > 2 had begun on trajectories with (p > (j), 
i.e. with a destabilised modulus, and thus access to G 
(see below) . We reiterate that we have only analysed the 
dS case in this example: such allowable cosmologies may 
undergo a different cycle of A/'ac before collapse in the 
AdS case. 



DISCUSSION 



A. Phenomenology of Fixed Points 



The analysis in Section (IV I showed the existence of 



many fixed points, with various degrees of stability. Dis- 
cussing the possible phenomenological implications of all 
of these would be a long and tedious process that we 
choose not to engage in. However, we will find it illumi- 
nating to discuss some properties of fixed point G, both 



long fluid dom. 




FIG. 12: Schematic summary of findings in {C, /3} space. C is 
the exponent in the modulus potential, and gives the coupling 
between axion and modulus. j3 is defined in Eq. ( 12 1 and 
represents a ratio of scales between the axion and modulus 
terms in the potential. They are the only two parameters 
that appear in the dynamical system, Eqs. (19 1. The vertical 



and horizontal lines for field oscillations are the conditions of 



The time scale for collapse out of B or G for similar trajectories 
with an AdS minimum would be similar. 



Eq. ([32| for ( ^ 0.3 dS vacuum (see Eqs. ([23}, (|24|), and they 
move in the directions shown for increasing (. The regions 
with A/'ac = 1 periods of accelerated expansion over all initial 
conditions only occurred for models from the FP-A simulation. 



by way of example, and since we will find them to be 
particularly interesting. 

The first thing to note about the fixed points is whether 
or not they occur in the bowl of the potential or on the 
plateau. We have that s/r = (j>/(f>, so that if s/r > 1 
the fixed point is on the plateau where the modulus is 
destabilised. For points G and F it is clear that they are 
on the plateau. Fixed point B has s = r = 0, so the 
ratio is undefined, and is technically at infinite modulus 
value. Hence in |49) temporary trapping was observed 
with large initial field values and a destabilised modulus, 
with exit from the fixed point occurring as the axion field 
value decayed. 

G is the only fixed point that allows for scaling solu- 
tions where both axion and modulus track the dominant 
fluid component. This was the tracking behaviour ob- 
served in |49j . This can be of particular use in allevi- 
ating fine tuning of axion initial misalignment angles in 
the following way. Heavy axions require flne tuning of 
their initial misalignment angle if they are not to "over- 
close" the universe by causing matter-radiation equality 
to occur at too high a redshift: they are outside the "an- 
thropic window" [HI [231 [Ml [HZ] ■ However, by allowing 
for tracking in the radiation era, the axion energy den- 
sity will scale as 1/a'*, instead of remaining a constant. 
Eventually oscillations will set in, since G is a saddle 
point, and the axion dark matter will scale as how- 
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ever this will begin from a lower energy density. The 
energy density is dumped into modulus kinetic energy of 
overdamped motion. This scaling will manifest as EDE, 
which we discuss in the next subsection. The difference 
to more standard tracking models is that the saddle point 
nature of G caused by the axion mass provides a natural 
mechanism for exit from tracking. Also in this model the 
tracking EDE field is not required to be the same as the 
field responsible for late time accelerated expansion, i.e. 
we have the additional pA- 



B. Early Dark Energy 

Temporary trapping in G during the radiation era is 
phenomenologically attractive because, as also pointed 
out in [49], it has the possibihty of leaving observable, 
and therefore constrainable, consequences as EDE. When 
is this situation possible? First, the modulus must be 
destabilised by axion initial misalignment, given by the 
bounds of Eqs. (131 and ( [T4| . The modulus will roll 
out towards x 



oo until Hubble friction stops it (if in 
additon C < \/6 then a temporary axion-modulus dom- 
ination in F will occur). Later, if the bound of Eq. (37) 
is satisfied, tracking will begin. The effects of this will 
further bound C. 

The axion-modulus EDE energy density contributes an 
amount He = ^<p + — '^Ib/C'^ during any period of 
tracking. EDE phenomenology places upper bounds on 
57e that translate simply to bounds on C: 



C > 29.- 



-1/2 



(42) 



During the radiation era the scaling EDE will behave 
as an extra effective relativistic species, AA'off , contribut- 
ing to the background expansion. The density contribu- 
tion can then be constrained by Big Bang Nucleosyn- 
thesis (BBN) and CMB bounds on iVgff. For example, 
the BBN constraints of [SB' allow for TVos = 3.85 ± 0.26, 
consistent with no change between BBN and the CMB. 
Taking the central value, parameterising the energy den- 
sity as [Ml 1100 ] , and assuming all the additional energy 
density to be in the form of EDE allows for 0,^ < 0.1: 



C > 6.2 BBN AiVoff' only. 



(43) 



This large value of fig would, however, be in conflict 
with the CMB (the agreement in [98] was for A'off only, 
and neutrinos behave differently in perturbations than 
EDE due to, for example, anisotropic stress). One of the 
main effects of the presence a sizeable Jig on the CMB 
is to change the location and amplitude of the acoustic 
peaks. The location of the first peak is related to the size 
of the sound horizon at decoupling which is given by 



da— 
da 



(44) 



Here, ^ = 3(1 -I- R) is the sound speed of the photon- 
baryon fluid and R{a) = |^ is the photon to baryon 



ratio. Using the Friedmann Equation today (subscript 
0) and at an epoch during the radiation era when the 
universe has evolved to point G yields: 



daV 
d7 j 



(45) 



We will assume that the presence of two background 
components, radiation, subscript 7, and matter (dark 
and baryonic), subscript m, does not change the result 



fie = 



C2 



during the radiation dominated era. Performing 



the integral in Eq. ( [44| (similarly to jlOll ) from the last 
scattering surface (Iss) to the epoch of matter radiation 
equality (eq) gives: 




(J2 



(46) 



where r^o is the standard sound horizon. The location 
off the flrst peak multipole is then: 



2tt 



C 



peak 



TsHo VC2 -4 



^0 I 



where the standard peak multipole is: 



^0 = 



27r 
rsoHo 



200. 



(47) 



(48) 



The qualitative behaviour is clear: for smaller C, i.e., 
for larger fig, the flrst peak occurs at a higher multi- 
pole. Ref. 1 102 performed simultaneous fits for fig and 
neutrino species along with other extended cosmologi- 
cal parameter sets, and found maximum values for f2e 
at the 95% confidence level of a few percent, with the 
absolute limit being dependent on priors about the DE 
(w > — 1: no crossing of the phantom divide) or the neu- 
trinos {Ncs > 3 from the standard model). The central 
value for N^s in these fits was around N^s ~ 3.6 ± 0.6. 
Taking the most generous upper limit of ilg < 0.042 
gives: 



C>9.8 CMB and AA^eff, [inij- 



(49) 



Currently, there is no detection of fig, but it will be 
possible to detect with current and future CMB experi- 
ments of Planck and CMBPol [SS]. Ref. [SS] reports, for 
a fiducial Planck central value of Jig = 0.03 and marginal- 
ising over their other extended DE parameters, an error 
of (Tfi^ = 0.003. A 3(7 measurement of fig translates to a 
bound: 



10-1 < C < 13.8 Planck forecast. 



(50) 



It is worth noting finally that the CMB and BBN 
bounds need not both apply, since the fields do not have 
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to have entered the scahng solution at any particular era, 
and can leave it. Of course there is also the caveat that 
these bounds only apply to the extent that motion in and 
near G is accurately described by the parameterisations 
used to derive them f |102| used a modified version of the 
parameterisation of |55j). and that approximately stable 
evolution in G can be maintained for long enough. 
In Fig. ( 12 ) we show a schematic for the phenomenol- 



ogy in different regions of {C, /3} parameter space that 
the results of this discussion and Section ( IV C I have led 
us to. 



The Assumption of Fixed fa, and Uplifting the 
Potential 



Throughout this work, as we have mentioned, we have 
assumed that fa can be taken fixed and that the modulus 
only effects the axion through exponentially scaling the 
mass. This had the simplifying property of providing a 
trivial metric on field space, with no change to the canon- 
ical kinetic terms. We can look at the validity of this 
assumption by computing what the change in fa would 
be along any particular trajectory. The assumption will 
be approximately valid if: 



|A/a|//a,o<l, 



(51) 



where fa.o is the point on the trajectory deemed to be 
"today" and the difference is calculated from the last 
relevant epoch. In the axiverse the scale of fa is fixed 
around lO^^ GeV by fixing the product 5* = Cx ~ 200. 
This does not appear in our dynamical systems analysis, 
since the scale of x only comes in via w, which the system 
does not depend on. 

As mentioned in Section [TTl] the change in fa will most 
likely be large for any trajectories that begin on the 
plateau of the potential and end in the bowl. This would 
require us to compute corrections in moving, for exam- 
ple, between fixed points B, G and the global minimum 
M. In the vicinity of the fixed points, checking that 
fa remains roughly fixed would require specifying uj and 
checking on a case-by-case basis. 

It is possible that trapping in G, or any fixed point 
with non-zero y, for an extended period of time could 
lead to large A fa/ fa- We can estimate this effect as 
follows. For a trapping of AA^ e-folds in G and setting 
Cxi = 200 to get the correct fa for the axiverse in the 
early universe, factors of C cancel and we have: 



Ax 



Afa ^ 37b 
/a, I 



200 



AiV = 0(10-2) AA^. 



(52) 



This will always be small for any scenarios of interest, 
since AA^ could only be large if G were driving infla- 
tion but we have seen that G itself cannot be acceler- 
ated and hence this is impossible, fa today will be only 
AA^% away from its initial value for small AA^, and there- 
fore predictions based on trapping in G in any particular 



epoch should be unaffected by our assumption of fixed 
fa- However, in predicting the fate of the universe, we 
emphasise again that G is unstable and moving into M 
in the future (or in the current epoch, as in [S'd') may 
entail large changes in fa- 
While the dynamical effect of changing the kinetic 
terms is hard to predict, it is simple to compute the 
change in the potential caused by identifying fa = 
The coupling term in the potential becomes: 



(53) 



This has one very interesting property: the emergence of 
a new, meta-stable (in the sense that it has a small barrier 
that can be tunnelled through, like the potentials of [7j) 
modulus minimum in the region of large (j). This meta- 
stable minimum can have positive cosmological constant, 
with no need for additional uplifting, i.e. with pA — 0. 
However, it is unstable in the axion direction, and could 
only usefully drive current accelerated expansion with an 
ultra-light axion of mass TOq < 10"^'^ eV. The emergence 
of the new minimum at large (j) can be traced to the extra 
term in d^U with opposite sign. 

For the new minimum to emerge one requires x < 2 /C 
which makes S ^ 0{1) and pushes fa — >■ Mpi. This leads 
to more fine tuning if this minimum is to provide late 
time acceleration since the small axion mass necessary 
for stability then needs to be put in by hand from the 
non-perturbative side, ruining the naturalness of the axi- 
verse scenario for light axions. We leave further study of 
the properties of this fa{x) scenario, particularly its pos- 
sibility of giving an alternative axion inflationary model, 
to a future work. 



VI. CONCLUSIONS 

In this paper we have studied a rich model of the dark 
sector, with many possible observational signatures as 
Dark Matter and Dark Energy, that extends and builds 
on well known work and tries to bring it into a broader 
theoretical context. Some of our findings are summarised 



in the schematic of Fig. ( 12 ) 



Axions and moduli are intimately linked to the prob- 
lem of the cosmological constant. Polchinski argued some 
time ago [S] , and indeed it has been known since the ear- 
liest days of string theory fSSJ I103| that the lightness and 
profusion of axions is a natural consequence of the the- 
ory, and is related to the anthropic demand for a small 
cosmological constant. Ultra-light axion fields with the 
hierarchy of masses generated by exponential dependence 
on the internal geometry of the compact space are obser- 
vationally relevant as a distinct form of dark matter. We 
have studied the cosmological evolution of axions when 
the energy scale of the potential is allowed to be dynam- 
ically controlled by a modulus of this geometry, instead 
of remaining fixed. If both the potential of the axion 
and the modulus arise from non-pertrubative physics at 
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similar energy scales, then we have shown that axion ini- 
tial misalignment can leave the modulus destabilised in 
the early universe and when the axion is allowed to be 
cosmologically light this can lead to significant evolution 
of the modulus throughout cosmic history. This evolu- 
tion allows for the possibility that the vacuum energy can 
change sign through the course of cosmic evolution, with 
today's quasi de-Sitter expansion being only temporary 
and the future evolution of the universe is a Big Crunch 
cosmic doomsday. We have demonstrated the topological 
change in the allowed phase space that such a possibility 
creates in two separate coordinate systems. 

By looking for fixed points of the dynamical system we 
have shown that a modulus destabilised by a large axion 
initial misalignment gives rise to a scaling solution where 
excess energy density in the axion field is redshifted away 
during the radiation era. One can always argue that a 
valid model must cut off the fine tuning on axion fields 
to values where there is a stable modulus, and we have 
shown the ratios of scales necessary for this. Alterna- 
tively, if destabilisation does occur and such a tracking 
behaviour ensues during cosmic evolution, then fine tun- 
ing on the axion dark matter is alleviated and in ad- 
dition observational limits on Early Dark Energy place 
constraints on the couplings of the model. 

We have not discussed the possibility of fitting this 
model to be cosmologically viable, which would require 
fixing Hq and w{z), among other things. Fits of this 
kind would allow comparison to current and projected 
constraints on w{z) and distance measurements (see, e.g. 
PUI 11041 1105| ). and in the context of this model trapped 
in B or G would allow predictions for future vacuum 
decay to dS in M 83^ or collapse [106l[l07]. It is worth 
noting, however, that even small uncertainties in the cur- 
vature, fife, can produce significant degeneracies and mis- 
estimations of w{z) from distance measurements alone 
|108| I109j , which highlights the need for more complete 
models, and use of more experimental estimators, when 
discussing non-standard models of Dark Energy. We have 
also not discussed perturbations, which would be neces- 
sary to compare this model properly to large scale struc- 
ture or CMB measurements. 

If the modulus in this model controlled a coupling of 
the standard model, then predictivity of any model build- 
ing will demand for it to be stabilised, and observational 
constraints will demand variations caused by axion evo- 
lution to be small, although potentially observable (see, 
e.g., [1101 llllp ^. Our analysis showed that this would 
lead to a tuning on axion initial misalignment, in addition 
to any related to dark matter density, if (pi < cj). If this 
bound is violated and the modulus is destabilised, any 



Many analyses of this kind, however, fail to account for the huge 
effect that variation in a would have on the standard model 
contribution to the vacuum energy through vacuum bubbles, at 
best greatly worsening fine tuning, and at worst ruining most 
anthropic explanations for the smallness of A. 



low energy constants that depend on it will be stabilised 
by Hubble friction and eventually scale according to the 
dynamics of a fixed point. We have shown that trap- 
ping in such a fixed point can reasonably maintain the 
axion decay constant, /q, and so may also be expected 
to naturally maintain any other constant with similar 
modular dependence. Anthropically, the meta-stability 
of this state of affairs is only as unnatural as a generic 
model allowing for a future Big Crunch. 

The future singularity allowed in the parameter space 
of this model changes the asymptotic structure of space- 
time and may be relevant to holographic models, or "Cos- 
mology/CFT" [TT2HTT5] . although a rolling rather than 
tunnelling to an AdS state in our model may trivialise 
any specific holographic mapping. Allowing for long lived 
unstable scalar potentials muddies the waters somewhat 
in the question of fine tuning in the landscape. The ax- 
iverse and supergravity (SU 11061 11161 1117) naturally al- 
low for scalar masses around Hq, but string quintessence 
models run up against many problems [118] , although for 
axions successful models do exist [85l Eg [T20]i°. How- 
ever, if the landscape favours instabilities |122j and as 
we have said they appear to be necessary feature in eter- 
nal inflation, it certainly seems pertinent to study their 
cosmology. Could it be that the seemingly unlikely situ- 
ation of many light axions pulling the moduli hither and 
thither in ultimately collapsing universes in fact opens up 
a whole new part of the "wasteland" , or that axion fric- 
tion can favour a large number of destabilised moduli and 
a natural route to non-trivial quintessence? What types 
of universe dominate the (admittedly controversially- 
defined) landscape volume: unstable, cosmological con- 
stant, cyclic or quintessence? We have also seen that a 
coupling of axions and moduli can allow for large varia- 
tions of the cosmological constant in the future, making 
multiple epochs of accelerated expansion possible during 
the matter dominated epoch in our past (observationally 
tightly constrained by ^^)^ or in the future evolution of 
the universe. Does this too effect our perception of fine 
tuning in relation to Dark Energy? 

In closing, we like to hope that the study of this model 
will motivate string theorists to further consider late time 
effects that the existence of ultra-light axions can have 
on diverse aspects of string cosmology, and demonstrate 
a new and rich model in the dark sector to cosmologists. 



During the final stages of preparation of this manuscript a very 
interesting model for natural, and indeed coupled, quintessence 
in string theory was proposed in 11211 . In particular, this involved 
a modulus controlling the size of a four-cycle, which in Type IIB 
theory can have a Ca axion associated to it. The mass of this 
axion will depend on the quintessence field, realising our model. 
In that work, the important constraints of fifth-force experiments 
and SUSY breaking are also addressed. 
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Appendix A: Stability Analysis 

In order to study the stability of the fixed points we 
expand about these points, setting X = Xc + (5X, with 
(5X the perturbations of the compact variables defined 
by Eqs. (18) considered as a column vector. To first 
order, the perturbations satisfy (5X' — W • (5X, where the 
matrix W contains the coefficients of the perturbation 
equations. The general solution for the evolution of these 
linear perturbations can be written as 



X = aiqie 



(Al) 



where the q„ are the eigenvectors associated with the 
eigenvalues A„ of the matrix W. Thus, the stability of 
the fixed points depends upon the nature of the eigen- 
values. The eigenvector of the corresponding eigenvalue 
determines the directions in phase space with which the 
eigenvalue is associated. We use the following classifica- 
tion [123] 

• (i) Stable node: A„ < for n = 1, 6 

• (ii) Unstable node: A„ > for n = 1, 6 

• (iii) Saddle point: A„ < 0, A™ > 0, (or A„ > 0, 

Xm < 0) for n,m — 1, ...,6 with n ^ m 

• (iv) Stable spiral: The determinant of the matrix 
W is negative and the real parts of the A„ are neg- 
ative 

For an expanding universe, a fixed point is an attrac- 
tor (stable) in the cases (i) and (iv), but it is not so in 
the cases (ii) and (iii) (unstable). We use the notation 
Xn{i, ■ ■ •} for an eigenvalue A. The subscript n labels 
the eigenvalue and the {«,•■■} denote the direction(s) 
in phase space with which this eigenvalue is associated, 
which may be determined by the n*^ eigenvector. 



The positive roots of the fixed points z^, r^, Sc, given 
in Table ( |Tl| lie in the expanding universe branch (-ff^), 
for which z,r, s,< > 0. The negative roots lie in the col- 
lapsing universe branch {H^), where z,r,s,t < 0. For a 
collapsing universe, a fixed point is stable if the eigenval- 
ues of W are positive. This is because the 'time' variable 
iV = In (a) of the autonomous system ( 19 1 becomes a 
decreasing function of time. 

The stability of the fixed points may be summarised 
as follows: 

• Point A This is the trivial solution corresponding to 
fluid domination where the kinetic and potential com- 
ponents of the axion and modulus fields plus A are 
negligible. It exists for all C, f3 and 7b. In H~^, the 
eigenvalues are 



Ai{x} = A2{y} = --(2-7b), 
XsM - \4{r} - A5{s} - A6{<} 



2 ' 



(A2) 



so this is a saddle point in the full phase space. Point 
A is unstable in the subspace of the {z,r,s,t} direc- 
tions and stable in the {x,y} subspace. In H~ this 
point remains unstable in the full phase space, since 
for realistic fluids 7b < 2. 

• Point B The p\ dominated asymptotic flxed point. 
In the eigenvalues are 



Xi{x} ^ X2{y} = -3, 

h{t} -37b , 

Xi{z} ^ X5{r} ^ Xe{s} ^ . 



(A3) 



Notice this point has three zero eigenvalues in the 
{z,r,s} subspace: We say this is a marginally stable 
solution in the sense that there is no instability growing 
exponentially, although it could be unstable to higher 
orders in the perturbation. To obtain the strict stabil- 
ity of this solution we would have to go beyond linear 
order, which we do not pursue as we know this point 
is ultimately unstable (see Section (IVB)). 



• Point C Corresponds to the axion and modulus ki- 
netic dominated solution. For this point to exist, 
— 1 ^ y !i 1- Rather than having an isolated fixed 
point, point C is formed of a continuous line of fixed 
points, known as an equilibrium manifold, which we 
call a critical line. This critical line is the unit circle 
+ 2/c — 1 ^^'l ^ symmetry of the autonomous sys- 
tem with Zc = rc — Sc = tc ^ 0. In both and 
the eigenvalues are 



Ai{z} = 3±V6Cy, 
A3W = 3, 



A5{s} = A6{0} = 3±^/^C2/, 



X2{x,y} = 6 - 37b ; 
X4{x,y} = 0, 



(A4) 
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where the ± in Ai and ^ corresponds to the ± of = 
±y. In General, if a nonhnear system has a critical 
line, the Jacobian matrix of the linearised system at 
a fixed point on the line has a zero eigenvalue with 
an associated eigenvector tangent to the critical line at 
the chosen point. Here, the zero eigenvalue A4 lies in a 



direction tangent to the 



Vc 



1 unit circle. This 



kind of non-linear system is a special subclass of the 

non-hyperbolic system (whose linearised system has 
one or more eigenvalues with zero real parts). 

The stability of a particular fixed point on the line can 
be determined by the non-zero eigenvalues, because 
near this fixed point there is essentially no dynamics 
along the critical line (i.e., along the direction of the 
eigenvector associated with the zero eigenvalue) , so the 
dynamics near this fixed point may be viewed in a re- 
duced phase space obtained by suppressing the zero 
eigenvalue direction. Then in this is point is un- 
stable in the full phase space since A3 > 0. 

In H~ , where stability corresponds to positive eigen- 
values, it might appear that stability is only guaran- 
teed for certain values of y along the fixed line. This 
is not so however, since both = +y and i/c = —y 
(with .Tc = ± ^/l — ) arc fixed point solutions in a 
collapsing universe and correspond to the upper and 
lower halves of the unit circle respectively. This guar- 
antees that every point along the critical line is stable. 
Hence, C is the asymptotic future of any model with 
an AdS vacuum. The particular fixed point along C 
that the system will finally evolve to will depend upon 
the initial conditions. 

We also note that As^e (which has a multiplicity of 2) 
is associated with only a single eigenvector, pointing in 
the s direction: the matrix of linearised perturbation 
coefficients is defective, i.e., it does not have a complete 
basis of eigenvectors, and is therefore not diagonaliz- 
able. As a result, the r direction is not represented 
and we so we write AejO}. This is almost certainly a 
consequence of the autonomous system having a higher 
dimensionality than required. 

Point D Modulus dominated fixed point, where Vb ^ 
Vd- This points exists for C < ^J~^. In the eigen- 
values are 



Ai{y,4=4(7'-37b, 
Xs{x} = 2C^-3, 
X,{r} = A5{s} = 



X2{y,z} = 2C^-3, 
A4{i} = 2C^ , 



To be stable along the Ai, A2,3 directions, which is the 

subspace {x,y,z}, requires C < and C < 

respectively. The {r,s,t} subspace is always unstable 
and hence point D is unstable in the full phase space. 
Point D remains unstable in the full phase space of 
H-. 

• Point E A scaling solution where then axion energy 
density vanishes and the modulus energy density scales 
with the dominant background fluid with Vb ^ Vd- 
This point exists for 7b < 2. In iJ+ the eigenvalues 
are: 



Ai,2{2/,4 = 



4C 



Chb - 2) 



±^(7b-2)(9C27b-67^-2C2) 



A3W = --(2-7b), 
3 

Asir-} = Xe{s} = -7b . 



A4{i} = : 



(A6) 



(A5) 



where the -|-root in Ai,2 is for Ai and the —root is for 
A2. For Ai,2 e ffi, C < yj^^^. Furthermore, if 
this condition is satisfied, A2 is negative V C, 7b and Ai 
is negative if C > \J~^^- When A 1,2 has a imaginary 
part, the real part is always negative. The {r,s,t} 
directions are always unstable and {a;} is stable for 
realistic 7b. We sec that point E is unstable in the 
full phase space and remains unstable in the full phase 
space of H~ . 



Point F Fixed point dominated by the modulus 

kinetic energy and potential energy of the axion- 
modulus coupling. This point exists for C < VQ. In 
the eigenvalues are 



Ai,2{2/,s} 



■7b- 



± - y'gC^ - 36(72 + 36C27b -|- 36 - 727b -|- 867^ - SC^jh - + SC^jh + 48(7 - 48Cjh , 
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AsW^^C-^-S, Xi{t} = lc\ (A7) 

X5{x,r} = 0, Xe{z} = -^C\ 

where the +root in A1.2 is for Ai and the —root is 
for A2. We see immediately that point F is unstable 
in the full phase space. In the presence of radiation, 
(7b = 4/3) there is a very finely tuned region of pa- 
rameter space open to C for which Ai and A2 may be 



real and negative. In the presence of dust, (7b = 1) the 
quantity under the square root simplifies to and the 
conditions for which Ai < and A2 < are C < \/3 
and C < \/6 respectively. In , point F is always 
unstable in the full phase space. 

• Point G A scaling solution where the axion and mod- 
ulus track the dominant background fluid. This point 
exists for 7b < 2. In the eigenvalues are: 



Ai,2{y,s} 



3 

4CP 



C2(7b-2)±yC(^ 



2473^7 



487^ 



- 24C72 + 8C27b + C37b - 2C3) 



^7b 



;7b: 



X6ix,r) = 0, 



(AS) 



where the -froot in Ai^2 is for Ai and the —root is for A2. 
For an expanding universe, there is a small region of 
C parameter space for both dust and radiation where 
Ai^2 are purely real and negative and a large region of 
C where A1.2 has a complex part with the real part 
negative. In these regions, {y, s} are stable directions. 
A5 is always positive however and so I is unstable in 
the full phase space of In , A3_4 are always 

negative and so I is unstable in the full phase space of 
H-. 

• Point I Modulus dominated fixed point, where Vd ^ 
Vb- Notice that this point mirrors point F and exists 
for C > a/6. In the eigenvalues are: 



A3,4{2^, s} 



2 12 ^ ^ 



6)(3C2 - 18 



8/32), 
(A9) 



Ai{y,r} = C2-37b: 



A2{y,r} 
Xe{z} 



-C^ , (AlO) 



where the -|-root in A3. 4 is for A3 and the —root is 
for A4. In H'^ regardless of the value of /3, the real 
part of A3 can only be positive. These constrains are 
due to the existence condition C > \/6. Hence, the 
{x, s} subspace is always unstable. The conditions for 
stability for Ai,2 in the expanding case are C < •\/37b 
and C < \/6, neither of which are ever satisfied if this 
point is to exist and so {j/, r} are unstable directions. 
We see that for both iJ+ and H~ , point I is unstable 
in the full phase space. 

• Point M The critical line in the {z, r, t} subspace 
corresponding to the global axion-modulus minimum. 



The emergence of this critical line is due to the fact 
that the autonomous system is of one dimension too 
many. With the minimally required number of au- 
tonomous system variables, the line would degener- 
ate to a unique point, given by Eq (30). M exists 



for all C , 7b , /3. It is important to point out that M 
is only a fixed point in the presence of a dS vacuum 
(C > 1/4) and corresponds to the asymptotic future. 
If the vacuum is AdS (C < 1 /4) the global minimum is 
not a fixed point, and the asymptotic future is cosmic 
doomsday in a Big Crunch. In the eigenvalues are: 



^±iv/9-24C2. 



Ai,2{y, 2,7-} 

A3{z, r,t} = -37b , 

A4,5{a:,s} = -^±^v/9- 8/3222 

Xelz^r^t} = 0. 



(All) 



where the -|-root in Ai_2 is for Ai and the —root is for 
A2. The -l-root in A4_5 is for A4 and the —root is for A5. 
Here, the eigenvector associated with Ag = points in 
a direction tangent to the fixed line and so the stability 
of a particular fixed point on the line can be determined 
by the non-zero eigenvalues, since near this fixed point 
there is essentially no dynamics along the fixed line 
(i.e., along the direction of the eigenvector associated 
with the zero eigenvalue) and the dynamics near this 
fixed point may be viewed in a reduced phase space 
obtained by suppressing the zero eigenvalue direction. 
Then, V C , 7b , /3 and for all points along the fixed line, 
the real parts of the eigenvalues are always negative 
and M is stable. There exist two bifurcation points 
along M: 



3 

8C2 



and z = 



9 

8/32 



(A12) 



Hence, for point M to be a stable node in the {y, z, r} 
subspace, z < •y/3/8C2, otherwise it is a stable spiral. 
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whilst for point M to be a stable node in the {x, s} 
subspace, z < -^/Q/S/?^, otherwise it is a stable spiral. 



Appendix B: 'Q'— Variable Autonomous System 



and the vacuum constraint C — -^pt- Suppressing the 
subscript i for brevity, the initial conditions are con- 
strained to lie on the three-dimensional manifold M.: 

^4 



= P-C- 



2 2 

r — z 



(CI) 



With the compact variables defined by Eqs. (26), the 
axion-modulus system may be written: 



Q 
v'q 



Q 
Q 



(I 

Q 



±Hc 



•cq - MsQy/zQtQ-^y , 



VQ 



C 



[2zl - ZQtQD^ + 4] , 



Q V2 



CVQ 



SQ + MxQy/zqt^ , 



with 



Q 



Q 



Hq 
2 



(Bl) 



+ 24 + 2y| 



/I 2 2 2 2 j.2\ 

7b(l ~ Xq -yg - Zq - Sq - Iq) 

T 



(B2) 



and 



HQ=±j3{l-ZQtQD^)^ 



(B3) 



where <& = 1/^/BpA- Here, (') 
is expanding, the sign of the square root in Eq. ( B3 ) 



^27- If the universe 



for Hq is positive. If the universe is contracting, the 
negative root should be chosen. 



Appendix C: Scanning the Initial Condition 
Manifold 



The autonomous system ( 19 1 has eight different pa- 



rameters which determine the subsequent motion of any 
given trajectory in phase space: six initial conditions, 
{xi,yi, Zi,ri, Si,ti} and two parameters, {C, /3}. To 
ensure that the entire space of initial conditions are 
sampled in an efficient around a region of interest, we 
use the method outlined in this appendix. 

We begin by acknowledging the Friedman constraint 

17b (initial) = l-{x^,+ yf + zl - rf + s\ + tf) , 



where p = \ — ^\, — x'^ — y^. For given values of x and y, 
the problem then reduces to varying two initial conditions 
evenly over A^, with the third constrained by Eq. (CI I. 



For example, if we were interested in scanning parameter 
space near to fixed point A, we could, with only some 
loss of generality, make the simplifying assumption that 
the axion and modulus fields begin frozen, a; = y = 0. 
For trajectories starting near to fixed point A, we have 
p = 1 — ilb ~ 0.01. We can then vary z and r, whilst 
still being free to independently vary {C, /?}. If on the 
other hand, we were interested in scanning parameter 
space near to fixed point G, since the fluid density and 
the value of y at the critical point depends on C, Jib = 



1 37b 



Vc 



I ^ , we are not free to independently vary 



C and {z ,r}. Each time C is changed, the equation for 
M also changes. This must be taken into account when 
using the method that is discussed below. 

The value of ( determines the topology of the initial 
condition manifold, which has non-constant curvature. 
If the vacuum is dS, C > |, the surface area of the man- 
ifold above some Smin is finite. For C ^ ji (AdS and 
Minkowski vacua) the manifold is not bounded above 
Smin and its surface area is infinite. This change in topol- 
ogy is most easily illustrated by considering the intersec- 
tion of the manifold A4 with the s = s^in plane. Solving 
Eq. ( |C1[ ) at s = Sniin for z, generates four roots, two of 
which are the ones applicable for this problem: 



f±ir) 



1 

71 



p- 



2pr^ — 2ps 



2 

mill 



4Cr4 



These two functions, f+{r) and f-{r) intersect at r = 
''max- Equating them and solving for r gives the physi- 
cally relevant root: 



'ip- 



i)(l + 2VC) 



(C2) 



V 4C - 1 

Until these functions intersect at 7'„iax we have /_ (r) < 



J+{r). Now, solving Eq. (Cll at s = Smin for r gives the 
two relevant roots: 



1/2 



which intersect at z = Zmax- Then, equating g+{z) and 
g-{z) and solving for z gives the physically relevant root: 



4C-1 



(C3) 
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Eqs. (C2) and ([C3j) illustrate the change in topology of 
M: for C = 1/4, r,„ax , Zmax oo and functions /±(r) 



and g± (z) and become parallel at large r and z and never 
meet; for ( < 1/4, r„iax , 2^max G C and the functions 
f±{r) and g±{z) diverge at large r and z. In both cases, 
the manifold never intersects the s — Smin plane. Only 
for C > 1/4 is the surface area of the manifold bounded 
For models where the vacuum is dS, if an 
is specified, the initial conditions on z and r may 



then be varied from some 



values given by Eqs. ( C2 1 and ( C3 ) 



mill ; ^min 



to the maximum 



Since M has non-constant curvature, it is not trivial 
to sample the manifold in a uniform way. To guaran- 
tee approximately uniform coverage, we project a grid 
of N squares, each of which have the same area. A, 
in a Cartesian coordinate system, onto Ai. Since Ai 
has non-constant curvature, the projected shapes will be 
four-sided polygons of unequal area. From the induced 
metric, gab, on M we can compute the surface area of 
each polygon. The components of gab read: 



511 

512 = 
522 = 



= 1 



1 

4^ 



4Cr3 



2r 



521 



4s2 



1 + ^ ^ - 2z 



The surface area Ai of each polygon i is then 



Ai 




'g dr dz , 



where 



9l2 



det [gab] = 511522 



+ z^p 



The limits of integration are the {r , z} boundaries of 
each square in the Cartesian coordinate system. If these 
limits lie outside of the manifold, the boundary functions, 
/±('")i where intersects Smin are used instead. 

The initial condition manifold is now divided into poly- 
gons of different area. We randomly 'throw' pairs of ini- 
tial condition coordinates {r^ , Zi} into each polygon on 
M . The number of initial condition pairs that are thrown 
is proportional to the area of each polygon: the larger the 
polygon area, the more initial condition pairs are thrown. 
Since the total number of points thrown onto the mani- 
fold is large, this statistical method ensures an approxi- 
mately even sampling of the initial condition manifold. 
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